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Abstract: 

We employ data taken by the JADE and OPAL experiments for an integrated QCD study 
in hadronic e"'"e~ annihilations at c.m.s. energies ranging from 35 GeV through 189 GeV. The 
study is based on jet-multiplicity related observables. The observables are obtained to high jet 
resolution scales with the JADE, Durham, Cambridge and cone jet finders, and compared with 
the predictions of various QCD and Monte Carlo models. The strong coupling strength, a^, is 
determined at each energy by fits of O(a^) calculations, as well as matched 0(q;^) and NLLA 
predictions, to the data. Matching schemes are compared, and the dependence of the results 
on the choice of the renormalization scale is investigated. The combination of the results using 
matched predictions gives 

as{Mzo) = 0.1187_o.ooi9- 

The strong coupling is also obtained, at lower precision, from O(a^) c.m.s. energy 

evolution of some of the observables. A qualitative comparison is made between the data and 
a recent MLLA prediction for mean jet multiplicities. 
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1 Introduction 

The renormalized strong coupling strength as of quantum chromo dynamics (QCD) is predicted 
to depend upon the momentum transfer of the interaction under study. It is desirable to perform 
tests of QCD at different values of this momentum scale. These tests are best performed using 
uniform experimental conditions since systematic uncertainties may thereby be reduced. Until 
recently, rather hmited sets of definite c.m.s. energies were available for QCD analyses in e+e~ 
collisions under consistent conditions. Since the start of the "LEP 2" program in 1995, the 
c.m.s. energy of the LEP collider at CERN has been increased in several steps from its original 
values close to 91.2 GeV, allowing QCD analyses over a wide range of high c.m.s. energies. 
However, the inclusion of measurements at lower c.m.s. energies is important as well, because 
QCD becomes more strongly dependent on the energy scale towards lower energies, and tests 
of the theory will be most significant here. 

This paper presents QCD tests using e"'"e~ annihilations into hadrons (so-called multi- 
hadronic events) from ^/s = 35 GeV to 189 GeV. Data recorded at the OPAL experiment 
at LEP are analyzed in combination with data from the JADE experiment at the PETRA col- 
hder at DESY, where e+e~ collisions were studied from 1978 to 1986 at lower c.m.s. energies. 
The OPAL and JADE detectors are similar in construction, and we have tried to keep the 
experimental procedure in both analyses as similar as possible. 

The observables used are exclusively based on the multiplicities of hadronic jets, defined 
using standard techniques. In the first part of the present work we present measurements of 
a large variety of such observables and compare them with several Monte Carlo predictions. 
The measurements are then employed to determine the strong coupling strength as and to 
test the QCD prediction for the momentum transfer dependence, i.e. the "running," of as- 
We compare the results from different types of matched O(a^) and NLLA predictions as well 
as pure 0{a1) predictions and study the dependence of the fit results on the renormalization 
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A 
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1986 
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1984/85 
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91.2 


1994 


34 


1508031 




130 


1995 


2.6 


144 


O 


136 


2.5 


140 


P 


130 


1997 


2.6 


179 


A 


136 


3.3 


167 


L 


161 


1996 


10.0 


281 




172 


1996 


10.42 


224 




183 


1997 


55.22 


1082 




189 


1998 


186.3 


3300 



Table 1: Overview of the data samples used in this analysis. The quoted numbers of selected 
events refer to the standard selection so described in Sect. ^. For the c.m.s. energies of 
the JADE experiment, the mean value {\fs)i averaged over the selected events, is given in 
parentheses. The small difference between the event number at \fs = 35 GeV and that given in 
is due to a somewhat more stringent definition of "long tracks" (see table |^) in this analysis. 

scale. The individual results are combined into a final value for as{Mzo). Finally, a recent 
MLLA calculation for the mean jet multiplicity is compared with our measurements without 
extracting a value for a^. 

2 The Experiments 

The e"'"e~ storage ring PETRA (see e.g. 0) was operated for physics measurements from 1978 
until 1986 at c.m.s. energies ranging from 12 GeV to 46.7 GeV. Extensive energy scans were 
made between these values. This analysis uses data samples recorded at energies of y/s = 35 
and 44 GeV. The precise ranges in c.m.s. energy and the luminosities are given in Table |^. 

The LEP e"'"e~ collider at CERN began operation in 1989 at ^/s ~ 91 GeV, i.e. around the 
mass of the boson. Since Fall 1995, the energy has been increased in steps from 130 GeV, 
136 GeV, 161 GeV, 172 GeV, 183 GeV to 189 GeV in 1998, the last year for which we include 
data in this study. The luminosities recorded by the OPAL experiment at each of these energies 
are also listed in Table ^ 

Descriptions of the JADE and OPAL detectors can be found in [|l], ^ and [Q, respectively. 
Apart from the dimensions, the detectors are very similar in their construction. Both are 
multi-purpose devices with a large solid angle coverage. Table ^ summarizes some detector 
parameters. The r, and z coordinates refer to a cylindrical coordinate system with the origin 
lying in the center of the detector and the z axis pointing along the incoming electron beam 
direction. 

The detector components primarily used in this analysis are the tracking systems and the 
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Parameter 


JADE 


OPAL 


Tl 1 TTl £in CI C 

j-^iiiitiiinK-iiin 


overall long 


th 




8 m 


12 m 


overall height 




7 m 


12 m 




dimension 


length 
outer radius 


2.4 m 
0.8 m 


4 m 
1.85 m 




transv. momentum 


At 


0.04 


0.02 




resolution (j{pt)/pt 


Bt 


0.018 


0.0015 




spatial 




r ~ <j) 


180 /^m (f 10 Atm) 


135 iiva 




resolution 
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1.6 cm 


4.5 — 6 cm 




double hit resolution 


7.5 mm (2 mm) 


2.5 mm 


Tracking 
system 


gas composition 
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OO . ( /O/ 0.0/0/ z . o /o 


oo/o/y.^/o/z.o/o 




gas pressure 


4 bar 


4 bar 




max. no. of hits 




48 
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reachable in 




0.83 • 47r 


0.73 • 47r 




at least 8 hits 
reachable in 


n 07 Am- 

u.y ( • 47r 


u.yo • 47r 




magnetic field 


0.48 T 


0.435 T 




energy 






0.015 


0.002 




resolution a{E)/E 


Be 


0.04 


0.063 




solid angle 


coverage 




90% 


98% 




angular resolution 


7 mrad 


2 mrad 






r 


adial extent 


1—1.4 m 


2.5—2.8 m 


Electromagnetic 






length 


3.6 m 


7 m 


calorimetry 


barrel 


polar angle covered 


32°-148° 


36°-144° 






radiation depth 


12.5Xo (15.7Xo) 


24.6Xo 








granularity 


8.5x10 cm^ 


10x10 cm^ 






outer radius 


0.9 m 


1.8 m 




endcap 


polar angle covered 


11°-327148°-169° 


11°-367144°-169° 






radiation depth 


9.6Xo 


22Xq 








granularity 


14x14 cm^ 


9x9 cm-^ 



Table 2: Parameters of the JADE and OPAL detectors are compiled from fT], ^, H. Wherever 
a second number is given in parentheses, it refers to a later period of detector operation. The 
significance of the parameters At, Bt-, Ac and Be is described in the text. 

electromagnetic calorimeters. The main parts of the tracking systems of both detectors are 
drift chambers built with a "jet chamber" geometry. The relative resolution of the transverse 
track momentum can be parameterized as (j{pt)/pt = a/A^ + {Bt ■ pt [GeV]y. The values of At 
and Bt for JADE and OPAL are given in Table 0. 

The electromagnetic calorimeters (ECAL) of JADE and OPAL are arrays of lead-glass blocks. 
The relative resolution of the electromagnetic cluster energy E is a{E)/E = A^+Bc/ ^/ E [GeV]. 
The parameters Ac and Be are given in Table ^ for both detectors. 

3 Data Samples and Multihadronic Event Selection 

In the JADE part of the analysis we use "detector level" information from data and Monte 
Carlo samples as described in [^], obtained from (measured or simulated) detector signals before 
any corrections are applied. We use data from the 1982 and 1986 runs at ^/s = 35 GeV and 
the 1985 runs at 44 GeV, where Monte Carlo samples including detector level information are 
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Reaction to be 
suppressed 


Cut variable 


JADE 


OPAL 


2-lcpton 
events 


nch 


> 3 long tracks 
and > 4 central tracks 


> ( 




Ileal 




> 7 




^shw 


> 3.0 GeV (barrel) 






or > 0.4 GeV (per endcap) 




2-photon 




> 0.5 


> 0.1 


events 


Pbal 


< 0.4 


< 0.6 




Praiss 


0.3 • 




other 


COS 9t 


< 0.8 


< 0.9 



Table 3: Comparison of the preselection cut values against two-lepton and two-photon events 
for the JADE and OPAL experiments, rich, ^^cai and Eghw denote, respectively, the number of 
charged tracks and calorimeter clusters and the total shower energy deposited in the calorime- 
ters. E^is is the total visible energy, Pmiss the total missing momentum and p^ai the momentum 
sum in z direction, normalized by Emg. Both experiments also apply cuts on the polar angle 
of the thrust axis, 9t- The terms "long" and "central tracks" are explained in Table ^ in 
Appendix |^. 



available. It was shown in previous reanalyses of JADE data 0] that measurements made in 
1986 and before could be reproduced using these data and Monte Carlo samples. 

At = M^o, we use the complete OPAL run of 1994 which is the largest available homo- 
geneous run without changes in the OPAL experimental set-up. 

From the first runs at higher c.m.s. energies in 1995 we combine the two samples at y/s = 130 
and 136 GeV, weighting them with their statistical errors. Runs in 1997 at the same two c.m.s. 
energies are used as well and subjected to the same treatment. The runs from 1995 and 1997 
are analyzed independently of each other since changes of the detector were made between the 
two dates. In both cases, results are quoted at y/s = 133 GeV. For each of the following energy 
steps we use the full event samples recorded by OPAL and quote results for the respective 
averaged c.m.s. energies. 

In order to reject lepton pairs and two-photon collisions, we apply standard selections from 
the two experiments. The selection used for the JADE part of the analysis is the same as 
described in and 0, and is closely related to that of OPAL, taking into account the difference 
in c.m.s. energies. For the OPAL analysis we apply the criteria given in 0. Table § contains 
a comparative list of the most important cuts used for the JADE and the OPAL analyses. We 
shall refer to this set of cuts as the "preselection" . 

At c.m.s. energies above M^o, photon radiation in the initial state becomes a significant 
source of background. In order to reject such "ISR events" we determine the total hadronic 
mass of an event following a procedure based on that described in [0 which takes possible 
multiple photon radiation into account. We require events to have y/s — y/s' < 10 GeV. For 
systematics studies we apply alternatively a combination of cuts on the visible energy and 



missing momentum of the event and on the energy of an isolated photon candidate |TT|. We 
shall refer to the former procedure as the "invariant mass" selection and to the latter as the 
"energy balance" selection. 

At a/s = 161 GeV and above, the production of (and later also pairs with hadronic 
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decays is another source of background. At these higher energies, we reject such reactions by 



dividing each event into hemispheres using the plane perpendicular to the thrust axis [|12[. We 



denote the heavier and lighter invariant mass of each hemisphere, normalized with the visible 



energy, by Mh and Ml [T^. Events with weak boson pairs decaying hadronically usually have 



larger hemisphere invariant masses. We apply the cut Mh/^+Ml < 0.35 and refer to this as the 
"jet mass selection". An alternative selection method has been used in previous OPAL analyses 
(e.g. IQ): the event is resolved into four jets using the Durham jet finder, and the 0{a1) QCD 
matrix element for a four-parton final state is calculated using the jet four-momenta pi...p4. 
The value of the matrix element is used as a cut variable. This method will be called the "event 
weight selection" and serves as a systematic check. The jet mass and event weight selections 
have very similar performance when applied to a sample of Monte Carlo events from PYTHIA 
[|T5| (for multihadronic events) and GRC4F (for all relevant four-fermion processes) which 
have been passed through the preselection and the invariant mass selection. 

Table ^ lists the numbers of events which we select at the individual c.m.s. energies using 
the standard selections. 



4 Measurement of Jet Fractions 

We present measurements of jet- multiplicity related quantities for various jet finders at c.m.s. 
energies of 35 through 189 GeV and compare them with predictions of several Monte Carlo 
models. The measured quantities are corrected for effects of limited detector resolution and 
acceptance, as well as for inefficiencies of the selection and ISR, i.e. they are presented at the 
"hadron level," which is understood to include all charged and neutral particles emerging after 
all intermediate particles with lifetimes below 3 ■ 10~^° s have decayed. 

4.1 General Analysis Procedure 

4.1.1 Reconstruction of Single Particles 

For both experiments, the measurements are based on charged tracks and electromagnetic 
calorimeter clusters. From these the four-vectors of single particles are reconstructed according 
to techniques which are quite similar in both the JADE and the OPAL part of the analysis. 
After imposing the quality criteria for charged tracks and electromagnetic calorimeter clusters 
described in Appendix ^ each accepted track is regarded as a charged particle having the 
measured three-momentum of the track and the mass of a charged pion. If a track can be 
linked to a particular ECAL cluster an estimate is made of the energy which a charged pion 
would deposit in the calorimeter; this amount is subtracted from the energy of the cluster. If 
the entire cluster energy is used up by such subtractions, the cluster is discarded; otherwise, 
the remaining energy is assumed to have been deposited by an additional neutral particle, and 
a zero-mass four-vector is constructed from this energy and the position of the cluster. 

The set of four-vectors from each selected event is then passed to the jet algorithms. At 
those c.m.s. energies where initial state photon radiation is large {^/s = 133 GeV and higher), 
the entire system of vectors is boosted into its own rest frame before the algorithms are applied. 

4.1.2 Correction for Experimental Effects 

At the c.m.s. energies of the JADE experiment (35 and 44 GeV), corrections rely on existing 
Monte Carlo samples with full detector simulation, taking into account the changes of the 
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detector with time. The detector level Monte Carlo samples available for this analysis were 
generated using the JETSET 6.3 program with the standard parameter set of the JADE 
collaboration |]18| and including photon radiation in the initial state. The Monte Carlo samples 
were shown in previous publications P, |^ to describe the data well. They were processed in the 
same way as the data, i.e. subjected to the same selection cuts and single particle reconstruction 
procedures. 

To obtain the corresponding hadron level information we have run the JETSET 6.3 program 
to generate events at both c.m.s. energies, using the same parameter settings, but without 
initial state photon radiation. By dividing the hadron level by the detector level prediction, 
binwise multiplicative correction factors were determined for all observables and then applied 
to the data. 

A similar binwise multiplicative correction procedure is applied to the OPAL data at ^/s = 91 
GeV. The correction factors were determined from two distinct Monte Carlo samples with 
detector simulation, generated by JETSET 7.4 and HERWIG 5.9 The parameter 



settings for both generators are described in |T^, ^ . Initial and final state photon radiation is 
included in both cases. 

For c.m.s. energies of 133 GeV and above, the JETSET Monte Carlo is replaced by PYTHIA 
5.7 [|T^, which has a more accurate modelling of initial state photon radiation. In addition, 
versions 5.8d (at y/s = 161 GeV) and 5.9 (otherwise) of HERWIG were used to generate 
multihadronic events. The correction procedure for ISR and detector effects applied at energies 
of 133 GeV and above was the same as at the lower c.m.s. energies. 

At y/s = 161 GeV and above, all relevant four-fermion final states were generated by the 



GRC4F generator ||T6| and subjected to detector simulation. The background predicted by 
GRC4F for each observable is subtracted at the detector level before the multiplicative correc- 
tions for residual ISR background and detector effects are applied. 



4.1.3 Determination of Systematic Errors 

To assess the size of systematic uncertainties inherent to the analysis procedure, the entire 
analysis was repeated with variations of the selection, of the correction Monte Carlo generators 
and of the detector components used. For each variation, the deviation of the final result from 
the standard measurement is taken as a systematic uncertainty. All systematic uncertainties 
are added quadratically to yield the total systematic error for every bin of each variable. 

The influence of the detector components used for the single particle reconstruction (tracking 
chambers and electromagnetic calorimeter) is estimated by repeating the analysis using only 
charged tracks. As a consistency check, where there are data sets from two different run periods 
{^/s = 35 GeV or 133 GeV), the influence of changes in the detectors carried out between the 
two dates has been investigated, but no effects were found. 

A systematic variation of the selection mechanism is done at all energies by tightening the 
cut on the thrust axis from |cos^r| < 0.9 to |cos0r| < 0.7. In addition, at ^/s = 35 GeV 
and 44 GeV, the cut variations described in are performed: The cut on the total missing 
three-momentum Pmiss is either tightened from 0.3 ■ ^/s to 0.25 ■ a/s or removed; the upper limit 
on the energy balance in beam direction, pbai, is either tightened from 0.4 to 0.3 or removed; the 
cut on the normalized visible energy Ems/\fs is varied from 0.5 to 0.55 and 0.45; a minimum 
of 7 rather than 3 "long tracks" (as defined in Appendix ^ is demanded. For the analysis 
at ^/s = 91 GeV, no further cross-checks aside from the variation of the thrust axis cut are 
done. At y/s = 133 GeV and above, the influence of the standard invariant mass selection is 
tested by replacing it by the energy balance selection, and at c.m.s. energies of 161 GeV and 
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higher, the analysis is repeated using the event weight selection rather than the standard jet 
mass selection. 

A variation of the correction Monte Carlo is not possible at JADE energies because no 
detector level samples aside from the ones used are currently available. At all other energies 
[s/s > 91 GeV), the analysis is repeated using HERWIG for the determination of the corrections 
instead of JETSET or PYTHIA. The effect of using the EXCALIBUR generator [||] rather 
than GRC4F for the background has also been investigated. The deviation from the main result 
induced by this change was found to be negligible throughout and is therefore not added as an 
additional uncertainty. 

Any systematic variation of the analysis procedure involving changes in the number of events 
entering the measurement of some numerical value (e.g. some bin) will necessarily generate a 
statistical deviation from the standard value. This effect becomes significant at c.m.s. energies 
with low statistics, i.e. at \/s = 133 GeV, 161 GeV and 172 GeV. To obtain a more realistic 
estimate of the systematic errors at these c.m.s. energies a number of subsamples are created 
from the respective multihadronic Monte Carlo sample, all of which contain on average the 
number of events corresponding to the integrated luminosity of the data (cf. [T^). Each of these 
subsamples is subjected to the same analysis procedure as the data. The standard deviation 
a of each separate contribution to the systematic error over all subsamples is determined for 
all measured values and regarded as the statistical component of the systematic error. If a is 
smaller than the corresponding systematic error contribution S measured from the data, the 
latter is reduced to y/S"^ — cr'^. Otherwise, the respective contribution to the systematic error 
is dropped completely. The procedure is performed separately for each contribution to the 
systematic error before they are added. The number of Monte Carlo subsamples used is 16 for 
the 1995 run at 133 GeV and 30 in all other cases. 

In order to further reduce fluctuations of the errors between measured values in adjacent 
bins due to low statistics, the overall systematic errors in all except the extreme bins of each 
quantity shown are averaged over each three adjacent values of that quantity. The systematic 
errors of the first and the last bin are subsequently set to the averaged errors of the three first 
and three last bins, respectively. This method of averaging systematic errors is performed at 
all c.m.s. energies. 



4.2 Description of the Measured Jet Fractions 



In the following, measurements are shown for three representative c.m.s. energy values: 35, 91 
and 189 GeV. All quantities are plotted versus the parameter(s) of the respective jet algorithm. 
Error bars represent the quadratic sum of systematic and statistical uncertainties. The latter 
include effects of limited statistics of the data and of the correction Monte Carlo samples. 
Listings of the numerical values, including those for the c.m.s. energies which are not shown in 
the figures, are to appear in the Durham data base [p^ . 

In each plot, the measurements are compared with hadron level predictions from the gener- 
ators PYTHIA 5.722, HERWIG 5.9, ARIADNE 4.08 H and CO JETS 6.23 HH, representing 
different kinds of parton shower and hadronization modelling. The parameter settings for ARI- 
ADNE and CO JETS are those from p4| and |2^. The predictions are based on samples of 
multihadronic events without initial state photon radiation, generated independently of those 
used for the corrections. In the ARIADNE samples, initial states generated by PYTHIA were 
subjected to the ARIADNE parton shower simulation, in which the hadronization modelling 
was taken over from JETSET 7.408. 
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4.2.1 n-Jet Fractions 



We analyze n-jet fractions, from tlie JADE |27|, Durham and Cambridge [2£] jet 
finders, wliicli reconstruct liadron jets based on different resolution parameters Ucut and different 
procedures to recombine unresolvable jets. Measurements for n = 2 through n = 5 and n > 6 
are shown versus ycut in Fig. 0, and in the upper and lower plots of Figs. |^ and |, respectively. 
For each variable, the measurements at the three c.m.s. energies are combined in one plot to 
make changes with energy more visible. The same ranges in abscissa and ordinate were chosen 
for all c.m.s. energies within one scheme. We also use a cone jet finder. 

At the highest c.m.s. energies, systematic effects due to hadronization are smallest. At 
^/s = 35 GeV and 44 GeV no HERWIG samples were available for the evaluation of the 
model uncertainties. Thus the apparent systematic errors are largest for the 91 GeV sample. 
The predictions of PYTHIA, HERWIG and ARIADNE are similar and lie mostly within the 
uncertainties of the measurements over the range in ?/c«t studied. It is impossible to choose a 
clear best generator from these three. Since measurements of jet fractions for neighboring y^ut 
values are partially based on the same events, the single bins of each jet fraction are strongly 
correlated. 

The predictions of COJETS are also in many cases in agreement with the data. However, 
especially at higher c.m.s. energies (i/i > 133 GeV), too many jets are predicted in the low y^ut 
regions. This is likely to be a consequence of the neglect of gluon coherence in the COJETS 
generator which may lead to an increased number of soft gluons between jets and therefore to 



high jet multiplicities. Coherence effects have already been observed by OPAL pO |. 

Scaling violations of both the Monte Carlo predictions and the measurements with y/s are 
visible for all three schemes in form of a reduction of the jet fractions for n > 3 (and a 
corresponding increase of R2) at higher ^/s, as is expected from the running of according 
to QCD. The resolution parameters ycut are defined in such a way that no scaling violations 
would be expected for jet rates if as did not run. 

Fig. I displays n-jet fractions obtained with a cone jet finder described in ||31|, which recon- 
structs jets within cones of half angle R having a minimum energy e. Because of the explicit 
energy cut inherent in the algorithm, results of the jet finding will depend strongly on the total 
energy E^g of the input particles. In order to remove this dependence, the cone algorithm is 
run with the scaled energy cut e' = e ■ Ems/\/s. In the upper part of the figure, jet fractions 
are plotted versus R for n < 2, n = 3 and n > 4 at a fixed minimum jet energy e. A value of 
e = 7 GeV is chosen for y/s = 91 GeV and above. At ^/s = 35 GeV and 44 GeV, where typical 
jet energies are considerably lower, e is set to 2.5 GeV and 3 GeV. The same jet fractions are 
shown against e in the lower part of the figures, where R is kept fixed at 0.7 rad. 

As in the case of the clustering algorithms, the predictions of all Monte Carlo programs 
except that of COJETS are similar for all c.m.s. energies. Again, COJETS deviates from them 
and from the data above ^/s = 91 GeV. 



4.2.2 yn Distributions and Differential n-Jet Fractions 

In the context of the JADE and Durham schemes, we shall denote by ?/„ the value of ycut at 
which a particular event switches from a n-jet to an n + 1-jet configuration^. Any one of the 
quantities ?/„ can be regarded as an event shape variable. We write differential distributions in 

^Strictly this definition is only reasonable provided that the yn fall off monotonically with n which is not 
always the case for the resolution measures considered here. 
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JADE Scheme 




^cut 



Figure 1: The fractions of events with 2, 3, 4, 5 and more than 5 jets as obtained by the JADE 
jet algorithm are plotted versus ycut at 35, 91 and 189 GeV. Data points are the measurements 
by the JADE and OPAL experiments. The error bars represent the total errors. Predictions 
by PYTHIA, HERWIG, ARIADNE and CO JETS are shown using hnes of different styles. 
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Durham Scheme 




Figure 2: In plots (a), (b), (c), (g) and (h), the fractions of events with 2, 3, 4, 5 and more than 5 jets as obtained by the Durham scheme 
are plotted versus Ucut at 35, 91 and 189 GcV. Diagrams (d) through (f) show the distributions in y„ (see Sect. 4.2.2) corresponding to 
(a), (b) and (c). The presentation is the same as in Fig. 1. Small horizontal ticks indicate the size of the statistical errors. 
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Figure 3: In plots (a), (b), (c), (g) and (h), jet fractions as obtained by the Cambridge scheme are plotted versus ycut at 35, 91 and 189 
GeV as in Fig. 2. Diagrams (d) through (f) show the differential jet fractions Dn (see section 4.2.2) corresponding to (a), (b) and (c). 
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Figure 4: The fractions of events with 2 or less, 3 and more than 3 jets as obtained by the cone algorithm are plotted in the upper 
row versus the cone half angle R at 35, 91 and 189 GeV, the jet energy cut-off e being fixed at 2.5 GeV for ^/s=35 GeV or 7 GeV 
otherwise. The lower plots show the corresponding fractions versus e with R kept fixed at 0.7 rad. Data and Monte Carlo predictions 
are presented in the same form as in figures 1 through 3. 



Vn as 

Y. ^ i^. (1) 

where a denotes the total hadronic cross section. In comparisons with theory, statistical corre- 
lations between the values at different parameter settings have to be taken into account. Such 
correlations are certainly present in n-jet fractions since the measurements at neighboring pa- 
rameter values contain partially the same events. In this respect, the differential distributions 
Yn are more convenient quantities since they contain each event only once. In the middle row 
of Fig. we present distributions Yn as obtained from the Durham scheme for n = 2 through 4. 

The Cambridge scheme employs a more complicated procedure for the jet reconstruction 
than the JADE and Durham schemes, involving two distinct resolution measures. There is 
therefore no counterpart of the quantities ?/„ in the Cambridge scheme which would allow the 
interpretation described above (see e.g. ^). We show instead, in the middle row of Fig. |^, 
the corresponding "differential n-jet fractions," 

A. ^ -P^, (2) 

with (T„ being the cross section for the production of n jets determined from explicit binwise 
differentiation of the jet fractions. The relation Dn = Yn — Yn-i, and in particular D2 = Y2, 
holds for conventional cluster algorithms like the JADE or Durham scheme, as long as the 
condition specified in footnote ^ is satisfied. 

Both for the Un distributions and the differential jet fractions, corrections and error calcu- 
lation were carried out independently of the measured jet fractions Rn- All plots show that 
the Monte Carlo predictions are almost indistinguishable and describe the data well except for 
CO JETS. None of the models is ruled out by the data on the basis of these quantities. There 
are scaling violations of the differential jet fractions as there were in the jet fractions. 



4.2.3 Mean Jet Multiplicities 

Another relevant quantity in the context of QCD tests is the mean jet multiplicity N defined 
by 

N=-y^nan, (3) 



a 

n 



for which various theoretical calculations exist. Measurements of versus ycut are presented 
for the Durham and Cambridge schemes in Fig. ^. The most complete theoretical predictions 
for this observable exist for these two schemes (cf. Appendix |B.1.2[ ). 

The decrease of the average number of jets at given ycut with rising c.m.s. energy, as predicted 
by the energy-dependence of the strong coupling, is clearly visible. Differences between the 
Durham and Cambridge results are predicted by all Monte Carlo programs and confirmed by 
the data. At and above ^/s = 133 GeV, the Durham scheme resolves systematically more jets 
than the Cambridge scheme. 

All Monte Carlo predictions except that of COJETS are almost identical and lie within the 
error bars of the measurements. The COJETS prediction overshoots the data at ^/s = 91 GeV 
and higher energies in regions of low ycut, i-e- in regions of large jet multiplicities, consistent 
with the explanation suggested in the previous section and with OPAL results from the analysis 



of charged particle based quantities |10, 11]. 
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Figure 5: The mean jet multiplicities as obtained by the Durham (a) and Cambridge schemes 
(b) are plotted versus i/cut at 35, 91 and 189 GeV. Data and Monte Carlo predictions are 
presented in the same form as in the previous figures. 



4.2.4 Mean values of ?/„ 

Finally, we consider the mean values of the ?/„ distributions given by 

iVn) = - [ Vn-^dyn (4) 

Measurements for n = 2 through 5 for the JADE and Durham schemes are plotted in Fig. |^, 
against the c.m.s. energy. Hadron level predictions of PYTHIA, HERWIG, ARIADNE and 
COJETS at each of the eight c.m.s. energies are also shown. To facilitate qualitative compar- 
isons, the respective predicted values at each c.m.s. energy are connected by spline functions 
to guide the eye. All generators are found in almost equal agreement with the data, except for 
COJETS. 

Higher moments of the observables yn have also been investigated. They have large uncer- 
tainties and are not presented. 



4.2.5 Summary 

The measurements of jet-multiplicity related observables show scaling violations. PYTHIA, 
HERWIG and ARIADNE describe the data well, but COJETS predicts too many jets at and 
above a/s = 91 GeV. 



5 Tests of Quantum Chromodynamics 

5.1 Procedure for ck^ Determinations 
5.1.1 General Description 

All jet-multiplicity based quantities considered in this analysis can be expressed as a power 
series in the strong coupling strength as, where the coefficients of the powers of as depend 
on the observable. The approximate methods used here for the calculation of such series. 
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Figure 6: The average of the observables i/n is plotted versus the c.m.s. energy ^/s for the JADE 
(left) and Durham schemes (right). Data points are the measurements by the JADE and OPAL 
experiments. The error bars represent the total errors. Small horizontal lines indicate the size 
of the statistical errors. The curves are spline functions connecting the Monte Carlo predictions 
at the single c.m.s. energies. 
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pure 0{a1) predictions and matched C(a^) and next-to-leading logarithmic predictions, are 
given in Appendix The coupling will be regarded as unknown quantity, and the tests of 
perturbative QCD will consist of fits of the predictions to the measurements over ranges in 
tjcut with as being fitted. These result in a determination of as at each c.m.s. energy. 

All fits are performed at the hadron level, using only the statistical errors for the calculation 
of x^- The error returned by the fit, i.e. the amount by which the fitted parameters may be 
varied without increasing by more than 1, defines the statistical errors of the fit results. 
Statistical correlations between bins in ycut are taken into account in the definition of x^- At 
c.m.s. energies where the statistics are sufficient (i.e. y/s=35 GeV, 44 GeV and 91 GeV) we 
determine the correlation matrix by subdividing the data sample, and the corresponding Monte 
Carlo sample used for the correction from detector to hadron level (see section 4.1.2 ), into 
independent subsamples and carrying out the entire analysis for each subsample. We choose 
AT = 30 at = 35 GeV, A^ = 15 at = 44 GeV and A^ = 80 at = 91 GeV. At the higher 
c.m.s. energies we determine the correlations from Monte Carlo subsamples. At ^/s = 133, 161 
and 172 GeV, 30 subsamples are used. At the two higher c.m.s. energies, we use 50 subsamples 
each. 

For the fits, we use only those observables for which there exist matched predictions of 0{al) 
and NLLA calculations, i.e. the 2-jet fractions and the mean jet multiplicities as measured using 
the Durham and Cambridge schemes (see Appendix 0). Fits are performed to the differential 
2-jet fraction D2 rather than the 2-jet fraction itself. Throughout the remainder of this chapter, 
the fitted variables shall be denoted D2, -D^, A^^ and A^*" where upper indices D and C stand 
for the "Durham" and "Cambridge" scheme. One aim of the analysis is a comparison between 
the different types of calculations. The functional expressions to be used in the fits are given 
for the various matching schemes in (p!5D , (|l^), ([T7| ) and ([l8|) of Appendix |B[ Furthermore, fits 
of the pure C(a^) predictions are carried out using (^ and (|^). The parameter to be varied is 
always as as appearing in the equations. The renormalization scale factor = /i/i/i is either 
kept fixed at = 1 or fitted simultaneously. Since the matched predictions are more complete 
than pure second-order calculations, one may expect an increased need for an adapted scale for 
the latter. 



5.1.2 Hadronization Effects 

In order to study the influence of hadronization on the quantities to be considered, we compare 
the predictions of PYTHIA 5.722, HERWIG 5.9 and ARIADNE 4.08 before and after the 
hadronization step. Each generator represents a different model for either hadronization or the 
partonic state. The generated event samples are the same as the ones used for the hadron 
level curves in Sect. [4.2| . The COJETS generator which uses independent fragmentation is not 
considered since it was seen to describe the data badly in certain kinematic regions. 

By the "parton level prediction," we understand quantities determined from the set of parti- 
cles emerging at the end of the parton shower generation, i.e. immediately before the hadroniza- 
tion step, including possible final state photon radiation. The curves in the lower partitions 
in Figs. ^ H and ^ are the ratios / of the parton level over the hadron level predictions for 
the differential 2-jet fractions and the mean jet multiplicities as a function of t/cut at all c.m.s. 
energies and for the three generators. The correction factors for PYTHIA, HERWIG and ARI- 
ADNE are shown, respectively, as solid, dashed and dotted lines. As is expected, hadronization 
corrections become notably smaller with rising c.m.s. energy. 

We transform the theoretical calculations to hadron level by dividing the calculation for each 
value of Ucut by the factors / shown in the figures. The result is then compared with the hadron 
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level measurement at the given Ucut- The PYTHIA generator will be used for the quoted central 
results. 



5.1.3 Determination of Fit Ranges 

The ranges of Ucut for the fits have to be chosen with care. A fundamental limitation is given 
by the kinematic region over which the respective theoretical predictions can be assumed to 
be valid. In particular, for pure 0{a1) fits, the choice of the fit range is limited to large Ucut 
regions. In the fits presented below, the fit range was adjusted separately for each type of 
calculation used. 

In order to be as independent as possible of assumptions on the hadronization process, we 
require hadronization corrections to be less than 10% in the case of the differential 2-jet rates 
and less than 5% in the case of the mean jet multiplicity and to be insensitive to changes 
in hadronization parameters and models. As can be seen from Figs. ^ to ^, these conditions 
have to be loosened at lower energies and for observables obtained with the Cambridge jet 
finder, in order to make as determinations possible at all. The condition of small hadronization 
corrections turns out to be the most stringent limitation of the fit ranges at low c.m.s. energies. 
The size of the corrections for experimental effects is also taken into account, but not considered 
as crucial because the simulation of the detectors is believed to be rather reliable and, in most 
cases, these corrections are clearly smaller than those for hadronization. 

As a further limitation we demand that the value of x^/d.o.f. obtained in the fit be not 
dominated by the contribution from a single bin at the boundary of the fit range. Another 
necessary check concerns the stability of the fit results under variations of the fit range. The 
results can only be regarded as reliable if they do not change significantly when the boundaries 
of the fit range are varied around the chosen values. 

Observing all conditions listed, we try to maximize the fit range in order to give the most 
significant possible results. 



5.1.4 Determination of Systematic Errors 

The determination of is repeated with variations in the details of the analysis procedure. 
For each variation, the absolute difference between the obtained value for and the central 
value is taken as the systematic error from the respective source. All contributions are added 
quadratically. 

At all energies, the same systematic variations of the selection procedure are performed as 
in the measurements of the jet quantities themselves (see Sect. |4.1.3D , and the fit is repeated 
using only charged tracks. At the c.m.s. energies of 35 GeV, 44 GeV, 91 GeV and 189 GeV, one 
observes in some cases a significant dependence of the resulting as value on the choice of the fit 
range. In order to estimate this uncertainty, both the upper and the lower boundaries of the fit 
ranges are varied by one bin in both directions, keeping the respective opposite boundary fixed. 
All contributions mentioned, added in quadrature, define the total experimental systematic 
error. 

Uncertainties from the hadronization models used for the transformation of the QCD predic- 
tion to hadron level are determined by varying the parameters and the Monte Carlo generators 
used. In particular, the parameter h appearing in the Lund fragmentation function, used for 
the hadronization of u, d and s quarks by the PYTHIA generator, and the width (jq of the 
transverse momentum distribution of the produced hadrons with respect to the parent parton 



are each varied in both directions within the uncertainties allowed by the OPAL tunes [21 1. At 
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y/s = 133 GeV and below, the parameter Qq defining the end of the parton shower cascade 
is also varied within its uncertainties around the central value of 1.9 GeV to 1.4 GeV and 2.4 
GeV, yielding, on the average, correspondingly smaller and higher numbers of partons at the 
end of the cascade. At the higher c.m.s. energies, these variations are replaced by a more 
radical change to Qo = 4 GeV, and the fits are repeated using HERWIG 5.9 and ARIADNE 
4.08 as described in Sect. 5.1.2| instead of PYTHIA. 



All theoretical predictions to be used make the assumption of zero quark masses. In an 
attempt to estimate the effect of this, we repeat the fit using a Monte Carlo sample containing 
only light primary quarks (u, d, s and c) for the transformation to hadron level. The quadratic 
sum of all contributions listed defines the overall hadronization uncertainty of the result. 

Finally, the uncertainty in the choice of the renormalization scale must be accounted for. In 
the fits where the scale factor is kept fixed at 1 for the central results, the fit is repeated 
with = 0.5 and = 2. Usually, the choice of a smaller scale will entail a decrease of the 
fitted as, while a larger scale will give larger values of as- The different deviations obtained 
from varying x^ in both directions are then added asymmetrically to the hadronization error 
to yield the total error from theory. In cases where both variations of x^ let as change in the 
same direction, the average of both deviations is taken and added as symmetric error. 

5.2 QCD Tests at Fixed c.m.s. Energies 
5.2.1 Fit Results 

Fits of the different matched calculations described in Appendix [B.1.2| , i.e. the InR matching, 
the R matching and their modified variants, as well as predictions have been performed 

for the measured observables D2, -Df , and A^*". In the case of the 0{al) calculations, fits 
both with a fixed QCD scale factor of 1 and with a variable scale were tried. The fitted 
predictions of the various matching schemes and the 0{al) calculations for ^/s = 35 GeV, 91 
GeV and 189 GeV are shown in the central parts of Figs. 0, ^ and ^ as smooth lines. Similar 
fits were done at all energies. At y/s = 35 GeV, no stable fits with a fitted QCD scale could be 
obtained for observables and . The small sections below each plot give the correction 
factors / for experimental and hadronization effects. In the additional small section above 
the plots, the differences between the fitted predictions and the data, denoted as S, is plotted, 
normalized to the measured value. The relative total and statistical errors on the data are 
shown along the 6 = line. 

The fitted values for as are summarized in Figs, [l^ and |lT] with their total errors (outer 
bars) and purely experimental errors (inner bars). The /d.o.f. is given for each result to the 
right of the plots. For all observables the same systematic pattern of the results from each 
calculation type used repeats at all c.m.s. energies. Furthermore, the theoretical uncertainties 
in the predictions resulting from the choice of the calculation type are at least as large as the 
experimental errors, in particular at y/s = 35 GeV, 91 GeV, 183 GeV and 189 GeV. Further 
refinements of the experimental procedure and increase of statistics will therefore not lead to 
significantly more precise results for as as long as theoretical uncertainties can not be reduced. 

Fits of 0{a'^) calculations with a fixed QCD scale of 1 require a limitation of the fit range to 
regions of large ycut in order to achieve acceptable values of x^/d-o.f. The resulting statistical 
errors are correspondingly large. As can be seen from Fig. |IU|, the results for as are generally 
high at lower c.m.s. energies, while, at ^/s > 133 GeV, they are similar to those obtained with 
the In R and modified i?-matchings. 

Allowing the scale factor x^ to vary allows the fit range to be extended to lower ycut- In 
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35 GeV 




Figure 7: The central parts of the plots display the results of fits of various theoretical pre- 
dictions (smooth curves) in comparison with hadron level measurements (data points) of the 
four observables D^, D^, and iV^ at = 35 GeV. The error bars on the data represent 
the total errors, the small horizontal lines the pure statistical errors which were used in the 
fits. Horizontal arrows indicate the respective fit ranges. The normalized difference between 
predictions and the data, S, is shown above each plot. The small sections below each plot 
indicate the correction factors / for experimental and hadronization effects. The hadronization 
corrections are shown as lines and shaded bands as described in the text. The experimental 
corrections are represented as open triangles and filled circles, respectively, for the 1982 and 
1986 runs. 



21 



91 GeV 




Figure 8: The fit results at y/s =91 GeV are shown in the same form as in Fig. |^. 
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Figure 9: The fit results at =189 GeV are shown in the same form as in Fig. 0. At this 
energy, the statistical fluctuations of the data points are too large to draw any meaningful 
conclusions from the normalized differences between the predictions and the measurement for 
D2 and D^- They are therefore not shown in the plot. 
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fact, a wider fit range is required in most cases to obtain stability of the fits. At the lower 
c.m.s. energies of 35 GeV and 44 GeV, the fit ranges have to be extended into regions where 
hadronization uncertainties become large, resulting in comparatively large theoretical errors of 
the results. 

In the case of the Durham scheme, the fitted as values obtained using 0{a1) predictions 
with a fitted scale are systematically smaller than the results from fits with a constant x^. In 
contrast, in the case of the Cambridge scheme they are somewhat larger. We conclude that a 
suitable choice oix^ is able to improve the agreement of the predictions with the measurements. 



a) Fit results for the 2-jet fractions 

The comparison of the quality of the fits (i.e. the /d.o.f.) does not indicate a clear preference 
for one of the matching schemes for either one of the variables. The normalized deviations 6 
also behave similarly for all matching schemes. Only in the case of the Cambridge scheme at 
\/s =91 GeV, they differ somewhat because rather distinct fit ranges had to be chosen for the 
different schemes. 

The results obtained with the In R and the modified In i?-matching with Umax = 1/3 turn out 
to be virtually identical up to at least the third decimal place of the extracted value for a^, which 
is in accordance with the findings of a previous OPAL analysis ||3^ . The fit results are therefore 
not shown separately, but only for the In i?-matching. The results of the modified /^-matching 
are also generally similar to those of the In i?-matching, but have slightly larger theoretical 
errors, /^-matching, however, leads to systematically low values of for both jet algorithms. 
In fact, the i?-matching results are in most cases the lowest of all calculation types under 



consideration. As is explained in Appendix |B.1.2| , the i?-matching may be expected to describe 



the data less well than the modified variant because the term G21L is not exponentiated. The 
observed behaviour seems to indicate that the inclusion of this term has more significance than 



the choice between Ini? and i?-matching. We decide to follow the practise of and use the 
In i?-matching for the final results. The difference plots 5 also indicate that the In R matching 
provides a good description of the data over a wider range in Ucut than both the R matching 
and 0{a1) calculations. 

The stability of the lni?-matching fits under variations of the fit range boundaries turns out 
to be generally good. 

b) Fit results for the mean jet multiplicities 

The fits of 0{a1) predictions for and A^*" with the QCD scale factor fixed at = 1 describe 
the data significantly less well in the regions of low Ucut than all the other fits, requiring a 
corresponding limitation of the fit ranges. Allowing to vary results in an improved agreement 
between data and predictions. Only at = 35 GeV, the 0{a1) calculations with x^ = 1 seem 
to follow the data well even in low Ucut regions. A closer inspection reveals, however, that the 
agreement at higher ?/c«t, i-e. within the selected fit range, is worse. 

The necessity of keeping hadronization uncertainties small limits the fit ranges rather strongly 



and leads to somewhat larger statistical errors than in the case of the 2-jet rates. Fig. [I^ shows 
that fits with a free QCD scale always result in values of as which are smaller than those from 
any other fits, while those with fixed to 1 yield usually the largest values. 

The agreement between data and theory, according to the normalized differences 5, is again 
rather similar for all matching schemes. All types of predictions display the property of over- 
shooting the data below some Ucut, i-e. predicting too many jets. 

For both jet algorithms, the fit results for as turn out to be less dependent on the choice of 
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Figure 10: Fitted values obtained using the different theoretical predictions described in 
the text for variable and D2 with the respective experimental (inner) and total (outer) 
errors. The x^/d-o.f. values of each fit are given in the rightmost column. 
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Figure 11: Fitted as values obtained using the different theoretical predictions described in the 
text for variable A^^ and A^*", presented as in Fig. |l^. At ^/s =35 GeV, no stable fits could be 
obtained using pure 0{al) calculations with a free QCD scale factor. 
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the matching scheme than is the case for the 2-jet rates, and the corresponding uncertainties 
are currently much smaller than the experimental errors. The modified In i?-matching of the 
Durham scheme contains the additional term LH21 in the exponent (see Appendix [B.1.2| ) and 



should therefore be a better description of the measurements than the In i?-matching. As can be 
seen from Fig. |Tl|, the results obtained by the modified In i?-matching, the i?-matching and the 
In i?-matching are always quite similar. Considering the small dependence on the calculation 
type mentioned above, one could expect to obtain rather precise determinations of from 
these variables if better hadronization models were available which would allow to extend the 
fit range towards lower i/cut- 

The mean jet multiplicities exhibit a stronger dependence on the choice of the boundaries 
of the fit ranges than the 2-jet rates, in particular at small c.m.s. energies. The fit range 
dependence (and therefore also the corresponding systematic error) is larger for the i?-matching 
and the modified In /^-matching than for In i?-matching. We conclude that the latter is a more 
appropriate description, valid over a wider range in i/cut, and decide to use it for the main 
quoted fit results. 

5.2.2 Investigation of the renormalization scale dependence 

The fits of the 0{a^) calculations show clearly that the data are better described by fixed 
order QCD predictions if small renormalization scales are used, with typically between 0.15 
and 0.2 in the case of the 2-jet fractions and between 0.03 and 0.1 in the case of the mean jet 
multiplicities. Similar results have been obtained in many other analyses involving event shape 



variables (e.g. 0, 0). The value of the best scale depends strongly on the observable. 
Although the fitted value of does not bear any physical significance, its deviation from unity 
indicates the importance of neglected higher order contributions in the fixed order calculation 
for the given observable. 

The inclusion of higher order terms in the matched predictions leads one to expect that the 
dependence of these predictions on the QCD scale factor will be reduced, as compared with 
the 0(0^) calculations. This has, in fact, been confirmed in |^, As was done in those 
pubhcations, we have studied the x^ dependence of the fit results by performing fits for at 
various fixed values of x^. Examples of the results for D2 and A^^ are shown in Fig. |l^ at 



•y/s = 35, 91 and 189 GeV. The solid lines demonstrate the behaviour of as as a function of 
Xfj, for both matched and 0{al) predictions. The fit ranges used for the matched predictions 
are the same as for the central values of as from the previous section. In the case of the 0(a^) 
predictions, the fit ranges from the simultaneous fits for ag and x^ are taken. The shape of the 
curves depends rather strongly on the fit range, which is why the circular and square markers 
are sometimes not situated on the corresponding curves. In general, an extension of the fit 
range causes the minima of the curves to become more pronounced, though the positions in 
x^ of the minima are found to be mostly unaffected by such changes. 

At 189 GeV statistics are low, and the fits of the C(a^) predictions tend to become unstable 
under variations of the QCD scale for larger x^. The corresponding curves are therefore not 
shown. Stability of the fits requires the limitation of the fit ranges to those used in the previous 
section for fits with fixed scales. Instabilities are further encountered at ^/s = 35 GeV in the 
case of the matched predictions of the 2-jet fractions at low x^, as well as for the two-parameter 
fits of the O(a^) calculations for the mean jet multiplicities. 



As can be seen in Fig. |I2|, the dependence of as on the scale in the vicinity of = 1 is 
generally much smaller with In i?-matching than with 0{al) calculations. Furthermore, fits 
of the matched predictions with a free QCD scale result in very shallow minima in and 
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Figure 12: Sensitivity of fit results to the scale x^. The solid and dashed hues are, respectively, 
the fitted and the corresponding /d.o.f. as a function of the QCD scale for observables 
D2 (left) and (right). Filled and empty squares denote and /d.o.f. for a fitted scale; 
circles have the same meaning for a fixed scale of 1. For at ^/s = 189 and at 35 GeV, 
results arc presented only for the In /?-matching scheme. In all other cases, fits of the pure 
0{al) predictions are also presented in the lower partitions of the plots. 



correspondingly large errors in the fitted parameters. Sometimes they do not converge at all. 
In some cases, convergence can be achieved by sufficiently extending the fit range towards 
low ycut- The results for are then usually found to be closer to 1. We conclude that the 
matched predictions indeed show little preference for any specific choice of the QCD scale, 
while the second-order calculations require rather definite values of to describe the data 
well. However, as can be seen in the figures, the Xfj_ dependence of the resulting as is still 
sizeable. 
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Figure 13: The combined results for as{Q) {Q = ^/s) from fits of the matched predictions 
plotted versus the c.m.s. energy. Arranged around the combined values are the contributing 
separate results for each of the four observables under consideration. All results are shown with 
their total errors (outer error bars) and the purely experimental (inner bars) component of the 
errors. A three-loop evolution of the current world average is overlaid as solid and dotted 
lines. 



5.2.3 Systematic Errors of the Main Result 

Based on the investigations carried out in the previous sections, we use the In i?-matching 
scheme for our final results for all observables. Tables ^ through [Tl| in Appendix y list the final 



results for for all variables and c.m.s. energies, along with the composition of their errors. 
We discuss the systematic errors for this scheme. 

At all c.m.s. energies, the use of only charged tracks rather than both tracks and electro- 
magnetic clusters has a large effect on the result. At JADE energies, variations of the selections 
generally contribute little to the systematic errors. At c.m.s. energies of 133 GeV and above, 
however, the influence of variations in the selections becomes sizeable. This is to be expected 
since the selections themselves play a more important role here than at the lower energies. The 
restriction of the cut on the thrust axis, which is applied at all c.m.s. energies, induces mostly 
small effects. 

The hadronization errors of the results are, in particular at lower c.m.s. energies, dominated 
by differences between the generators (i.e. HERWIG or ARIADNE). This could be anticipated 
considering the size of the hadronization corrections shown in Figs. ^ to ^ A similar behaviour 
is seen for the dependence of the results on which is largest at low c.m.s. energies. The 
influence of the omission of b quarks in the hadronization model and of varying the termination 
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0.0030 
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±0.0014 
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0.0025 


± 


0.0016 
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0.1188± 0.0045 ± 
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±0.0030 
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Table 4: Combined results from fits of matched predictions, along with the statistical, experi- 
mental systematic, hadronization and QCD scale errors (left column). The values in the right 
column are results of a three- loop evolution of each result to Q = M^o. 



point of the parton shower is sizeable and relatively independent of the c.m.s. energy. Variations 
of the hadronization parameters b and ag have very small effects throughout. 



5.2.4 Combination of the Fit Results 

We combine the results obtained using the four observables, separately at each c.m.s. energy by 
forming the mean value, taking into account the covariance matrix corresponding to the total 
errors of the individual values. The statistical correlations are determined from the data/Monte 
Carlo subsamples of Sect. The systematic uncertainties are treated as fully correlated, 

aside from the variations of and those of the bin boundaries which are considered as uncor- 
related and added to the diagonal of the matrix only. To this end, the scale uncertainties are 
made symmetric by taking the mean value of the two separate deviations. In order to break 
up the total error into its single components, the mean is determined using, respectively, the 
covariance matrix corresponding to only statistical, only experimental, and only experimental 
and hadronization uncertainties. Fig. |l^ gives an overview of the combined results, as well as 
the results from each single observable with their experimental and total errors^. The combined 
results are given in Table ^ The values listed in the right column of the table are three-loop 
evolutions to y/s = M^o. The results are compatible within their errors. 



As can be seen in the Fig. 13, the separate results are distributed in roughly the same pattern 



at each c.m.s. energy. The solid and dotted lines represent a three- loop evolution of the current 



world average of ^^(M^o) = 0.119 ±0.004 ||3^ with its total error. No significant discrepancies 
between the curve and any of the fit results can be noted. Our values agree with the world 
average within their experimental errors. Our results using the 2-jet fraction of the Durham 
scheme are also in good agreement with previous determinations using resummed calculations 
at these energies (see, e.g., §, |T0|, |n| g, |3|] ; 



The difference Aa^ = ctsiQ) — as{Mzo) is plotted in Fig. |T4[ Again, the inner error bars in 



■^At ^/s = 172 GeV the mean value lies above all contributing individual values. Here, relatively large positive 
correlations and the fact that the highest individual value (from N^) has the smallest error have the combined 
effect that the lowest is obtained when all individual values lie on the same, rather than on opposite sides of 
the mean. 



30 



QCD with a (M„)=0.119 



J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 



_L 



-0.02 
-0.03 

20 40 60 80 100 120 140 160 180 200 

Q [GeV] 

Figure 14: The difference Aa^ = as(Q) — as{Mzo) versus Q = a/s from ct^ fits of matched 
predictions. The definition of the error bars is as in Fig. |13[ The sohd hne represents the 
current world average. 



the plot represent the total experimental errors while the outer bars include also hadronization 
and QCD scale uncertainties. The statistical and experimental systematic errors are obtained 
from those quoted in Table § assuming that correlations between different energies can be 
neglected. In order to take correlations of the theoretical errors into account, each difference 
Aag was calculated separately under each variation of the hadronization model and of the QCD 
scale. As usual, the difference between the value of Aa^ obtained after each such variation and 
its central value was counted as a contribution to the systematic error. The contributions to 
the hadronization error were combined in quadrature. 

Finally, in order to obtain a single result for a^, we perform a fit of the three-loop running 
of as to the values of Fig. |13|, the fitted variable being ^^(M^o). We assume that correlations 
between the total experimental errors can be neglected and use these for the determination 
of x^- The resulting fit error is taken as the total experimental error of the final result. The 
statistical component is assessed by repeating the fit using only statistical errors. The total 
hadronization and scale errors are then calculated as for Aa^, i.e. by repeating the fits under 
each variation. The final result of 0.1187laooi9 is given with its individual error contributions 
and the value of x^/d-o.f. in Table |[ 



5.3 as{Mzo) from the Energy Evolution of Jet Observables 

The availability of measurements performed over a wide range of c.m.s. energy with very similar 
experimental conditions suggests investigations of the c.m.s. energy dependence of jet-related 
observables rather than analyses at fixed energies. In Sect. |, the energy evolution of jet 
fractions was already qualitatively observed, and the extracted as values shown in Fig. |T^ were 
seen to be in accordance with the QCD prediction of the running of the strong coupling. 

Similar analyses have been done based on data from many e^e~ collision experiments, present 
and past. For an early test of the running of as, the three-jet fraction was used (e.g. This 
is, in lowest order, simply proportional to a^. The JADE scheme was employed for these 
analyses because it has the smallest and least energy- dependent hadronization corrections of 
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Table 5: The result for as{Mzo) obtained from a fit of the three- loop running of the strong 
coupling to the combined values determined using matched QCD calculations at single c.m.s. 
energies (Table with the individual error contributions. For each systematic variation of 
Monte Carlo parameters and of the QCD scale, the deviation with respect to the central value 
is given with a sign indicating the direction of the deviation. 



all clustering schemes at low c.m.s. energies. More recent measurements have been mainly 
concerned with the measurements of the energy evolution of event shape moments and of the 
mean jet multiplicity (e.g. H, [11, |3^, 



In this analysis, the energy evolution of the three-jet fraction at a fixed value of ycut and 
that of the mean values of the observable y2 are compared with the predictions of QCD for 
various clustering algorithms. 

5.3.1 Analysis Procedure 

The comparisons between measurement and predictions are done by making fits to the re- 
spective observables over the entire range of c.m.s. energies considered in our analysis, taking 
as{Mzo) as the fitted variable. The expressions for the 0{as) and 0{a1) coefficients in the pre- 
diction are taken from |^ for the JADE scheme and from ||2^ for the Durham and Cambridge 



schemes. The expression to be fitted to the mean value of ?/2 is given by (|^) of Appendix 
The coefficients A and B were obtained from runs of the EVENT2 program pTj 

In the case of the fits to the three-jet fractions, a specific value of Ucut has to be selected 
at which the jet fraction is given. Here, we select a value in a region where the sensitivity to 
changes in the c.m.s. energy is large (see Figs. |1|, ^ and 0) and hadronization corrections are 



small. In the case of the JADE scheme, the value of ycut = 0.08 from is used. For the 
Durham and Cambridge schemes, we take ycut = 0.01. 

For the moments of some event shape variables, there exist non-perturbative corrections 
(also called "power corrections") to the perturbative calculations It was found in |^ that 



corrections of order l/Q as well as 1/Q to the mean values of y2 as obtained with the Durham 
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Figure 15: Mean values of y2 as obtained with the JADE and Durham schemes at parton level 
versus the c.m.s. energy ^/s. The solid lines represent the fitted 0{al) predictions at a fixed 
QCD scale of 1. 



scheme can be neglected at the c.m.s. energies considered in this analysis. Hadronization 

corrections for the thrcc-jct fractions at the chosen value of ycut are found to be small (between 
2% and 11%) for the JADE and Durham schemes, but large (about 25%) for the Cambridge 
scheme at low c.m.s. energies. We perform the fits for all observables at parton level. The 
transformation of the measurements from hadron to parton level is done separately at each 
c.m.s. energy by means of multiplicative factors which are obtained from the same Monte 
Carlo samples as were used for the fits. The factors predicted by the PYTHIA 5.722 Monte 
Carlo at standard OPAL tune are used to correct the results. 

The systematic errors can not be treated by repeating the fits with purely statistical errors 
for each systematic variation, because the systematic variations applied are essentially different 
at each c.m.s. energy. We therefore use the total experimental errors of each measurement for 
the calculation of the ^-nd take the resulting fit error as the overall experimental error of 
as- Its statistical component is determined by repeating the fit using only statistical errors. 
Hadronization uncertainties are determined by repeating the fits with variations of the Monte 
Carlo predictions for the correction factors from hadron to parton level. We apply only the 
variations leading to the predominant error contributions in the fits of the previous section, i.e. 
changes to HERWIG and ARIADNE and the removal of the b quark. The absolute differences 
from the central result are again added in quadrature to get the total hadronization error. 

Fits are carried out both with a constant QCD scale factor, fixed at unity, and with taken 
as an additional fit parameter. In the first case, systematic variations of to 0.5 and 2 are 
performed and added asymmetrically to the hadronization error. 
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Figure 16: The 3-jet fraction as obtained with the JADE, Durham and Cambridge schemes at 
parton level versus the c.m.s. energy ^i. The solid lines represent the fitted O^a^) predictions 
at a fixed QCD scale of 1, the dashed lines are the results obtained with a fitted scale. 



5.3.2 Results from the Energy Evolution Fits 



Fig. |T5| shows the measurements of (1/2) for the JADE and the Durham scheme at parton 
level with the result of the fit at fixed x^. The inner error bars denote the size of the purely 
experimental errors, the outer bars include also hadronization errors. The obtained values for 
as(M^o) from the fits shown are listed in Table 12 of Appendix with the deviations induced 



by each systematic variation, as before. The total errors are of the same order of magnitude as 
from the fits at separate c.m.s. energies. Fits with a free QCD scale lead to errors of almost 
100% in the scale for both observables. We therefore present the results obtained at a fixed 
scale of 1, performing the usual variations to estimate the error from the scale uncertainty. The 
somewhat small values of x^/d.o.f. indicate correlations in the systematic uncertainties. 

The results of the fits to the three-jet fractions are summarized in Table |13| of Appendix |C| 
and are displayed in Fig. Unlike the previous fits to the mean values of ?/2; the fits with a 
free QCD scale lead to smaller errors in x^. The corresponding prediction is added as a dashed 
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as{Mz») 


0.1181 


Stat, error 


±0.0006 


Total cxp. 


±0.0020 


udsc only 


+0.0015 


HERWIG 


-0.0035 


ARIADNE 


-0.0011 


Total hadronization 


±0.0039 


= 0.5 


-0.0035 


— 2 


+0.0048 


Total error 


+0.0066 
-0.0056 



Table 6: The weighted mean of as{Mzo) as determined from the energy evolution of observables 
{112)1 (2/2^)5 -^3; -Rf and i?^, with the single error contributions. For each systematic variation 
of Monte Carlo parameters and of the QCD scale, the deviation with respect to the central 
value is given with a sign indicating the direction of the deviation. 



line in the plots. The optimized scales again turn out to be significantly smaller than = 1. 
The fit results for with a fitted scale are systematically smaller than with a fixed scale of 
1. In order to be able to estimate the error induced by scale uncertainties, we quote central 
results for x^ = 1 and vary the scale to 0.5 and 2. 

We combine the separate results with = 1 from the five observables into one value for 
as{Mzo)- The combined result is calculated by taking the average of the remaining five single 
values, weighted with their respective total errors. We assume total experimental errors to be 
largely unaffected by correlations and determine these also by taking the weighted mean of the 
single values. The error contributions from each hadronization and QCD scale variation are 
determined in the same way separately for each variation, resulting in the overall hadroniza- 
tion and asymmetric scale error of the combined value. Table ^ shows the combined result 
O.llSl^aoose with its single error components. 



5.4 MLLA Prediction for the Mean Jet Multiplicity 

Lastly, we compare our measurements of the mean jet multiplicity of the Durham scheme with 
a recent hadron level prediction [H3|. The calculation was carried out in the framework of the 
modified leading logarithmic approximation (MLLA) in the form of a cascade of successive 
parton branchings which is continued until the relative transverse momentum between the two 
partons emerging from the branching falls below a cut-off Qq. Assuming "local parton-hadron 
duality" to be valid and setting Qo to a value of typical hadron masses, the prediction may be 
compared with hadron level measurements. The lower portions of each part of Fig. |1^ show 
the hadron level measurement of at all c.m.s. energies with their errors as in Fig. |^ and 
the MLLA prediction |^ as a solid line. The cut-off Qq is fixed by In {Qq/ Kqcd) = 0.015 with 



^QCD = 500 MeV This value was obtained from fits to measured hadron multiplicities 
over a c.m.s. energy range from 1.5 GeV to 91 GeV. The smaller captions above each plot show 
the normalized difference S between measurements and predictions, defined as in Sect. f).2.1[ A 
first observation is that the MLLA curves rise more rapidly than the data below some ycut as a 
result of the singularity in the strong coupling as the transverse parton momentum comes close 
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Figure 17: The mean jet multiplicity as obtained with the Durham scheme plotted versus 



parameter y^ut at all energies. The MLLA prediction according to for each energy is the 
solid line and the same prediction shifted by Ql/s in ycut as described in the text is the dashed 
line. 
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to Aqcd- The difference seen at low values of i/cut becomes smaller for higher c.m.s. energies. 

The MLLA prediction assumes quarks to be massless. In order to artificially introduce 
a hadron mass in the calculation, it is suggested in |^ that one may compare the parton 
transverse momenta k± with the transverse energies of the hadrons given by E± = \/kj_ + Q^. 
In the application of the predictions to the Durham jet multiplicities, this amounts to comparing 
the measurement at some ycut with the prediction at y^ut + Qo/^^ i-^- a relative shift by a fixed 
amount along the ycut axis. The dashed lines in the plots represent the prediction after this 
shift. At ^/s = 35 GeV and 44 GeV, the predictions somewhat undershoot the data after the 
shift, showing that such a crude method for introducing a hadron mass into the calculations 
works less well at such low energies. 

At all c.m.s. energies, the data fall below the MLLA predictions in the region of medium 
ycut- This behaviour is also observed in [BH] for = 91 GeV and is there attributed to the 



omission of higher loops in the definition of as 



6 Summary 

We have performed QCD related measurements at c.m.s. energies of 35, 44, 91.2, 133, 161, 172, 
183 and 189 GeV using data from the JADE and OPAL experiments which are very similar in 
their components used in the analysis. The same measuring techniques were applied in the two 
experiments, including details of the selections and the method of reconstructing single particle 
four-momenta. The results of the measurements display the same systematic behaviour over 
the entire c.m.s. energy range, suggesting that the desired homogeneity of the analysis is indeed 
realized. 

Hadron level measurements of n-jet fractions, differential jet fractions, distributions in the 
variables ?/„ and mean jet multiplicities have been presented up to large jet multiplicities and 
down to very low values of Pcut using the JADE, Durham and Cambridge jet finders. Mea- 
surements of jet fractions as obtained with a cone algorithm, varying both of its parameters R 
and e, have also been performed. The mean values of the ?/„ distributions have been measured. 
The numerical values will appear in the Durham data base. The measured values were com- 
pared qualitatively with the predictions of four Monte Carlo generators, PYTHIA, HERWIG, 
ARIADNE and COJETS, representing the major currently available models for parton shower 
evolution and hadronization. All generators except for COJETS were found to be in agreement 
with the data. COJETS was seen to predict too many jets, in particular in regions of high 
jet multiplicities (high jet resolution). The discrepancy rises with c.m.s. energy and becomes 
significant at ^/s > 91 GeV. It can be explained by the omission of gluon coherence effects in 
the generator which will lead to an excess of soft gluons. Qualitatively, all observables based 
on the clustering schemes display a clear scaling violation with c.m.s. energy, as expected from 
QCD. 

The 2-jet fractions and mean jet multiplicities as obtained with the Durham and Cambridge 
schemes were used for quantitative tests of QCD. Matched C(a^) and NLLA predictions for 
these observables were fitted to the data at each separate c.m.s. energy over an appropriate 
range of ycut, the fitted parameter being the strong coupling at the respective energy. The 
InR and the i?-matching schemes, as well as their "modified" variants were used in the fits. 
In addition, O(a^) predictions were fitted, where the QCD scale factor was taken as an 
additional free parameter or kept fixed at unity. The fit quality in terms of was found to 
be reasonable, confirming the general validity of QCD within the errors of the results. None 
of the different calculation types could be clearly disqualified on the basis of the values. 
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To obtain an acceptable x^, the fits of the 0{a'^) calculations at = 1 required a limitation 
of the fit range to regions of large i/cut, which is in accordance with the expectation that the 
omission of higher orders of becomes noticeable in regions of high jet multiplicities (low 
Ucut)- The range of validity of the 0{al) calculations can be extended towards lower ycut by 
choosing an appropriate scale x^. In agreement with other analyses, the optimized values for 

turn out to be significantly lower than 1. At all c.m.s. energies, the resulting ag values 
display the same relative shifts as a function of the type of prediction chosen. The results using 
O(a^) calculations and = 1 tend to be systematically larger and the corresponding results 
with a fitted scale tend to be smaller than those from the matched predictions. In the case of 
the 2-jet fractions, the differences between the results using simple and modified In i?-matching 
are neghgible. The results of the modified iti-matching are also generally comparable with 
those of the In i?-matching, while the /2-matching, which is known to be less complete, leads 
to systematically small values of a^. In the case of the mean jet multiplicities, the differences 
between the matching schemes are less significant. Generally, we observe that at c.m.s. energies 
with sufficiently high statistics (35 GeV, 91 GeV, 183 GeV and 189 GeV) the precision of the 
obtained ag is currently limited by the theoretical uncertainties. 

The dependence of the results on the QCD scale has been investigated. Simultaneous fits 
of as and x^, using matched predictions and the same fit ranges as with a fixed scale either 
do not converge at all or result in values for x^ around unity affected by very large fit errors. 
The as results, however, still depend on the scale, which leads to sizeable contributions to the 
theoretical errors. 

Combinations of the In i?-matching results from the four observables, separately at each 
c.m.s. energy, are found to be in agreement with a three-loop QCD evolution of the current 
world average for ^^(M^o) of 0.119 ± 0.004. A fit of the three-loop running expression for ag 
to the combined results over all c.m.s. energies returned a final value of 

X „TTe„+0.0034 
as{Mzo) = 0.1187_o.ooi9 

with a x^/d.o.f. of 7.85/7. 

We have also carried out fits to the energy evolution of the mean values (1/2) as obtained 
with the JADE and Durham schemes and of the 3-jet fraction of the JADE, Durham and 
Cambridge schemes, each evaluated at a fixed tjcut- In all cases, ^^(M^o) was fitted for, and 
predictions of 0(q!s) were used with a fixed QCD scale of unity. The results found are again 
in agreement with the world average with reasonable x^/d-o.f., except for the 3-jet fraction of 
the Cambridge scheme where a somewhat high value of about x^/d.o.f. = 31/7 is obtained. 
The differences between the results from different observables turn out to be larger than the 
overall experimental error. Attempts to fit the QCD scale simultaneously resulted in very large 
fit errors in the case of the quantities (1/2), precluding any reasonable fixing of the scale. In the 
case of the 3-jet fractions, optimized QCD scales were again found smaller than unity and led 
to systematically smaller results for ^^(M^o). The weighted mean value of the ag results of the 
five observables at = 1 is as{Mzo) = O.llSl^aoosg. 

Finally, we have tested a hadron level MLLA prediction for the mean jet multiplicity in 
the Durham scheme, separately at each c.m.s. energy. A qualitative comparison with the 
data showed that the prediction overshoots the data in regions of medium ycut which may 
be attributed to the fact that as is included only in one-loop accuracy. Another significant 
deviation is seen towards very low ycut, where the prediction begins to rise significantly faster 
than the data due to the singularity in the ag running. This deviation is largest at low c.m.s. 
energies. Better agreement between the data and the prediction in these regions can be achieved 
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if the prediction is shifted in Ucut by an amount corresponding to typical hadron masses. 

In summary, this analysis presents a unique investigation of the running of in a large 
c.m.s. energy range of 35 through 189 GeV based on a consistent treatment of the data and 
employing up-to-date theoretical predictions. The numerical value of ag obtained from our 
study of jet rates and jet multiplicities is found to be in agreement with the world average, 
which has been obtained from a large variety of observables and processes. It is of comparable 
precision. 



Acknowledgement s : 

We particularly wish to thank the SL Division for the efficient operation of the LEP accelerator 

at all energies and for their continuing close cooperation with our experimental group. We 

thank our colleagues from CEA, DAPNIA/SPP, CE-Saclay for their efforts over the years on 

the time-of-flight and trigger systems which we continue to use. In addition to the support staff 

at our own institutions we are pleased to acknowledge the 

Department of Energy, USA, 

National Science Foundation, USA, 

Particle Physics and Astronomy Research Council, UK, 

Natural Sciences and Engineering Research Council, Canada, 

Israel Science Foundation, administered by the Israel Academy of Science and Humanities, 

Minerva Gesellschaft, 

Benoziyo Center for High Energy Physics, 

Japanese Ministry of Education, Science and Culture (the Monbusho) and a grant under the 
Monbusho International Science Research Program, 
Japanese Society for the Promotion of Science (JSPS), 
German Israeli Bi-national Science Foundation (GIF), 

Bundesministerium fiir Bildung, Wissenschaft, Forschung und Technologic, Germany, 
National Research Council of Canada, 
Research Corporation, USA, 

Hungarian Foundation for Scientific Research, OTKA T-029328, T023793 and OTKA F-023259. 



A Quality Criteria for Charged Tracks and Electromag- 
netic Clusters 

The quality cuts defining acceptable charged tracks are listed for the JADE and OPAL experi- 
ments in Table |^. For the JADE part of the analysis, two different types of tracks, "long" and 
"central" tracks, are defined. In the case of OPAL the criteria used in the preselection differ 
from those applied in the actual analysis. Uhus denotes the number of hits in the respective 
central jet chamber, and Uexp the number of hits which is to be expected taking into account 
track direction and detector geometry. For the OPAL preselection hits in the entire central 
detector were counted. The symbol p denotes the reconstructed three-momentum of a track, 
its projection onto the xy plane and 9 its polar angle. The beam energy dependent upper 
limit Pmax{Ebeam) is taken to be (1 + 6a/0.022 + (0.0015Efeeam)^)-E6eam- The upper limits on p 
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JADE 




OPAL 




Variable 


"long 


"central 


preselection 


analysis 


analysis 




tracks" 


tracks" 




(tracks only) 


(tracks & ECAL) 


^hits 


> 25 


- 


> 20 


> 40 


> 20 








— 




> 0.5 


p [GeVJ 




> 0.1 




< 90 


Pmax\-^beam) 


P± [GeV] 


> 0.5 




> 0.05 


> 0.15 


> 0.15 


\zq\ [cm] 




< 35 


< 40 


< 25 


< 30 


do [cm] 




< 3 


< 2 


< 2 


< 2.5 


cos 6' 






< 0.995 




< 0.9622 


i?i [cm] 






< 60 







Table 7: Quality criteria for charged tracks at JADE and OPAL 



Detector 


Variable 


JADE 


OPAL 


region 






preselection 


analysis 


barrel 


Edust [GeV] 


> 0.15 


> 0.1 


> 0.1 








> 1 


> 1 


endcap 


Edust [GeV] 


> 0.15 


> 0.2 


> 0.25 








> 2 


> 2 



Table 8: Quality criteria for electromagnetic calorimeter clusters at JADE and OPAL 

are motivated by the occurrence of left-right ambiguities for tracks running close to and parallel 
to a wire plane. In some cases, this situation leads to the reconstruction of very straight tracks 
with extremely high momenta. The expression for Pmax{Ebeam) is the beam energy, augmented 
by six times the track momentum resolution, zq is the z coordinate of the point of closest 
approach (p.c.a.) of the fitted track helix to the origin of the coordinate system, and do and 
Ri are the transverse distances from, respectively, the p.c.a. and the first hit to the origin. 
The quality criteria for electromagnetic calorimeter clusters are given in Table ^ In the case of 
OPAL, the criteria used in the preselection differ again from those applied in the actual anal- 
ysis. Edust and rihi are, respectively, cluster energy and number of lead-glass blocks contained 
in the cluster. 

B Theoretical Calculations for Jet-Multiplicity Related 
Observables 

B.l Second-Order Approximations 

Predictions up to fixed second order in ag are available for the jet fractions and mean jet 
multiplicities of all jet finders presented in Sect. ^ The perturbative power series for the 2-jet 
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fractions and the mean jet multiplicity can be written in terms of = a^/ {2tt) as 



Ft. 



cut ) 



1 + AR{ycut)as + BR{ycut)al 



(5) 



and 



cut I 



2 + AN{ycut)as + l3N{ycut)al (6) 

where indices "i?" and 'W on the coefficient functions indicate the 2-jet fraction and the 
mean jet multiplicity, respectively. The coefficient functions A and B may be calculated by 
integration of the O(a^) matrix elements. We employ second-order predictions for the 2-jet 
fractions and mean jet multiplicities as obtained with the Durham and Cambridge schemes for 
tests of perturbative QCD. For the Durham scheme, we use values for the respective coefficients 



obtained for a wide range in ycut |0] from numerical matrix element integration using the 
program EVENT2 In the case of the Cambridge scheme, the calculations from valid 
for 0.001 < ycut < 0.2, are used. 

predictions for the mean values of ?/„ can be obtained for n 



respectively, of next-to-leading and leading order: 



(1/2) 



Aas + Bat 



2 and n = 3, being, 



(7) 



(^3)«("?) = Ca', (8) 

Very recently, calculations of four-jet observables in next-to-leading fixed order (i.e. 0{al)) 
approximation have been presented E^. We do not use these for QCD tests, but we present 



measurements of the relevant quantities in Sect. 



B.1.1 Next-To-Leading Logarithmic Approximations 

The truncation of the perturbative series after a fixed order is meaningful as long as the omitted 
terms can be assumed to be small. For the jet-related observables under consideration, as well 
as for cumulative cross sections of event shape variables, the coefficients can themselves be 
written as series in L = ln(l/?/cut), which becomes large in the region of small ycut- For some 
observables, reasonable predictions of cross-sections may nevertheless be obtained in these 
kinematic regions of low ycut, because they allow the "resummation" of the largest logarithmic 
components of the coefficients to all orders in as Both the differential 2-jet fraction and 
the mean jet multiplicities as obtained with the Durham and Cambridge schemes belong to this 
group of observables, which is why we chose to use them to test QCD. 

In the case of the differential 2-jet rates, the feasibility of resummation is closely connected 
to the fact that the complete prediction can be written in the "exponentiated" form 

R2{ycut) = C{as)expG{as,L) + D{as,L) (9) 

where 

00 

c(«,) = i + ^a«: (10) 

n=l 

and 

00 n+l 
n=l m=l 

= Lgi{asL) + g2{asL) + asQsiasL) + alg4{asL) + ... (11) 
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and the remainder function D{as-, L) is assumed to vanish for ycut ~^ 0. Exponentiation imphes 
that terms with m > n + 1 are absent in the sum (pA]) . As can be seen from ([lID , each of 
the functions Qi is defined as an infinite power series in a^L. The two functions, gi and g2, 
contributing, respectively, the largest and second-largest logarithmic parts of the coefficients, 
have been calculated for the 2-jet fractions of the Durham and Cambridge schemes |^^, yielding 
a prediction in "next-to-leading logarithmic" approximation (NLLA). Up to this order, the 
predictions are identical for both schemes [^. The coefficient Ci can be obtained from the 
0{as) matrix element and C2 from integration of the O(ttf) matrix element and comparison 



of the result with (^, e.g. by fitting the latter to the first 

An improvement of the NLLA prediction suggested in concerns the subleading loga- 



rithmic term G2i<y'^L, which is included in the second-order prediction and may therefore be 
determined from a comparison of the NLLA with the 0{al) prediction as was done for C2. A 



result is again quoted in The overall NLLA prediction is then 



Ri'-'-^iycut) = (1 + Cia2 + C^al) exp + g2{asL)] . (12) 

NLLA predictions have also been derived for the mean jet multiplicities of the Durham scheme 
50| . Here, however, the complete prediction does not exponentiate as in @. One therefore 



has to include powers of L above n + 1 and obtains 

00 2n 



n=l m=0 

„2n , r-1 



= 2 + h^{Lal) + L-^h2{Lal) + L-'h^{Lal) + (13) 

where the "leading" and "next-to-leading" logarithmic terms are now those with m = 2n and 
m = 2?T, — 1, respectively. Each function hi is an infinite series in La^. The complete NLLA 
predictions are given in [Q. 

As in the case of the 2-jet fraction, one may additionally obtain the coefficients of subleading 
logarithmic terms of 0{as) by comparing them with the fixed order prediction. Here, the terms 
in question are -^22-^^ and H21L. In fact, from expansions of the calculations in an analytic 



expression for H22 can be derived Determinations of both H22 and if 21 have been performed 



using fits to Bn HJ, 0- Including these two terms, the overall NLLA prediction takes the form 
iV^''''^(2/c«t) = 2 + /ii(La2) + L-^h2{Lal) + H22L'' + H21L. (14) 
Again, the prediction turns out identical for the Cambridge scheme. 



B.1.2 Matched Predictions 

Various techniques have been devised to combine fixed-order with logarithmic predictions in 
order to obtain a prediction which is valid over a wide range of Ucut- In the most obvious of 
such combination schemes, generally called the "i?-matching" scheme, the exponential function 
in (0) is expanded in its power series and the coefficients ci and C2 of the as and a1 terms are 
read off. The matched prediction R2{ycut) is then formed by subtracting the 0{as) and 0{a1) 
terms from R2^^^ and adding the second-order prediction, yielding 

R2{ycut) = R^'-'-^iycut) - c^a, - C2al + R^'^"''\ycut). (15) 

The procedure requires the explicit knowledge of Gn, G12, G22 and G'23 whose values can be 
obtained by expanding function G{asi L) (with the known functions gi and (72) in a power series 
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and reading off the appropriate coefficients. We shall follow the convention of [^] and use the 
expression "i?- matching" if the G21 term is left out in i?^^^"^ (and consequently also in C2). If 
the term is kept, we shall speak of "modified i?-matching" . Because of the additionally included 
subleading term, one may expect the modified variant to be the more precise prediction. 

In the corresponding matching procedure for the mean jet multiplicity the 0{as) and 0{al) 
coefficients which appear in both types of predictions can be read off directly from (|13D resulting 
in the matched prediction 



N^^^^{ycut) - diois - C2<^1 + AN{ycut)ois + BN{ycut)ai 



cut I 



(16) 



It does not make any difference in this case whether the subleading terms -^22-^^ and H21L are 
included or not, because they are subtracted out in any case. The required coefficients Hn, 
H12 

m 



H23 and H24: have been calculated analytically and are given in equation (8) of |50|. 
An alternative way of matching the two types of calculations, called "In i?-matching" in 



applies a similar procedure as before, but now to the logarithms of the predictions. To this end, 
the logarithms of the two predictions are expanded in their power series and terms of 0{a^) as 
well as non-logarithmic terms are omitted. The 0{as) and O(a^) terms are then subtracted 
from the resulting NLLA expression, and the two logarithms are added to yield the matched 
prediction 



lni?2(2/, 



cut I 



Lgi{asL) + g2{asL) 

-{GuL + Gi2L')as - {G22L' + G23L')al 

1 



(17) 



BR{y, 



cut I 



:^Riy, 



cut ) 



The subleading term G2ia'^L in the NLLA part of the prediction can be ignored since it is 
removed in the process, being implicitly contained in Br. 

It has been pointed out in that neither the i?-matching nor the lni?-matching scheme 
ensures that the resulting prediction for R2 tends to 1 in the kinematic hmit of ycut = ymax 
where all events are resolved in two jets. In the case of the In i?-matching, this drawback can 
simply be cured by replacing L with L' = ln(l/?/cut — ^/ymax + !)• As in |Q, we refer to this 
variant as "modified In i?-matching" . 

The mean jet multiplicity may be subjected to a procedure analogous to the In i?-matching 
where the logarithm is taken of — 1 rather than A^ itself. The expansion of ln(A^-^^^'^ — 1) 
makes it possible to read off the common 0{as) and 0(a;^) coefficients di and d2. Combining 
both predictions and subtracting double terms leads then to the matched result 



HN{y,^t) - 1] 



ln[Ar^^^^(y,„,) 



1] 



diGs - d2al + AN{ycut)0is + 



BN{y. 



cut) 



Here, it does make a difference whether the subleading terms -^22-^^ and H21L are left out. 
We shall speak of "In i?-matching" if only the analytic expression for H22 is included, and of a 
"modified In i?-matching" if the fitted values for both H22 and H21 are used. 

R and In /^-matching differ generally in their assumptions on the unknown subleading and 
non-logarithmic terms of 0{al). In the case of the i?-matching, these terms are ignored in the 
predictions for R2 and A^, while, in the case of the In /^-matching, they are assumed to vanish 
in the logarithms. 
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C Fit Results for ag 



The tables on the subsequent pages summarize the values obtained for «<, in Sects. |5.1| and 
p.3| and the detailed composition of the systematic errors. For each systematic variation, the 
deviation with respect to the central result is given with a sign indicating the direction of the 
deviation. 













a^(35 GeV) 


0.1482 


0.1438 


0.1510 


0.1500 


Fit range [log^g (yc«t)] 


-1.55 — -0.65 


-1.55 — -0.65 


-1.4 — -0.5 


-1.55 — -0.65 


xVd.o.f. 


5.25/5 


7.66/5 


0.27/6 


0.66/6 


Statistical error 


±0.0011 


±0.0013 


±0.0039 


±0.0029 


Experimental syst. 


±0.0031 


±0.0036 


±0.0057 


±0.0043 


Total hadronization 


±0.0073 


±0.0047 


±0.0046 


±0.0063 


Xfj, = 0.5 


-0.0003 


+0.0007 


-0.0038 


-0.0019 


x^ = 2 


+0.0055 


+0.0044 


+0.0068 


+0.0054 


Total error 


+0.0097 
-0.008 


+0.0066 
-0.0066 


+0.0107 
-0.0091 


+0.0098 
-0.0084 













a^(44 GeV) 


0.1294 


0.1257 


0.1353 


0.1313 


Fit range [log^o (yc«t)] 


-1.55 — -0.65 


-1.85 — -0.95 


-1.7 — -0.5 


-1.85 — -0.65 


xVd.o.f. 


4.21/5 


12.36/5 


3.56/8 


4.58/8 


Statistical error 


±0.0013 


±0.0016 


±0.0028 


±0.0017 


Experimental syst. 


±0.0066 


±0.0053 


±0.0074 


±0.0062 


Total hadronization 


±0.0028 


±0.0067 


±0.0036 


±0.0056 


Xfj, = 0.5 


-0.0001 


+0.0004 


-0.0024 


-0.0004 


x^ 2 


+0.0036 


+0.0030 


+0.0047 


+0.0031 


Total error 


+0.0081 
-0.0073 


+0.0088 
-0.0088 


+0.0099 
-0.0090 


+0.0091 
-0.0085 













a^(91.14 GeV) 


0.1214 


0.1131 


0.1218 


0.1163 


Fit range [log^oiy cut)] 


-2.125 — -1.375 


-1.5 — -0.625 


-2.25 — -0.625 


-1.625 — -0.5 


xVd.o.f. 


8.52/5 


29.44/6 


11.69/13 


1.59/9 


Statistical error 


±0.0002 


±0.0002 


±0.0003 


±0.0008 


Experimental syst. 


±0.0038 


±0.0051 


±0.0037 


±0.0056 


Total hadronization 


±0.0032 


±0.0021 


±0.0027 


±0.0022 


x^ = 0.5 


-0.0006 


+0.0007 


-0.0013 


-0.0009 


x^ — 2 


+0.0035 


+0.0019 


+0.0032 


+0.0026 


Total error 


+0.0061 
-0.0051 


+0.0056 
-0.0056 


+0.0056 
-0.0047 


+0.0066 
-0.0061 



Table 9: Results for a.^ from fits of the In /^-matching predictions for the four observables 
, D^, and A^^ at = 35, 44 and 91 GeV 
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7V^' 


a, (133 GeV) 


0.1066 


0.1047 


0.1153 


0.1126 


Fit i-Migc [logx„((y,;„()] 


-2.5 -U.3 


-2.25 -0.3 


-2.3 -0.5 


-2.23 -0.5 


xVd-o.f. 


5.69/7 


6.1/6 


5.63/8 


5.63/7 


Statistical error 


±0.0021 


±0.0035 


±0.0028 


±0.0042 


Experimental syst. 


±0.0076 


±0.0057 


±0.0078 


±0.0062 


Total hadronization 


±0.0024 


±0.0011 


±0.0021 


±0.0017 


•r,, = 0..J 


<().()()()1 


+0.()00() 


-0.0010 


-U.0003 


x^ = 2 


+0.0021 


+0.0015 


+0.0027 


+0.0020 


Total error 


+0.0085 


+0.0068 


+0.0090 


+0.0079 


-0.0082 


-0.0068 


-0.0086 


-0.0077 













as{l6l GcV) 


0.1097 


0.0991 


0.1106 


0.1090 


Fit range [logio 


-2.75 0.75 


-2.5 0.75 


-2.75 0.5 


-2.75 0.5 


\-/d.o.L 


<S.().S/7 


i:-!.42/(i 


3.9(i/9 


4.!)/9 


Statistical error 


±0.0038 


±0.0039 


±0.0028 


±0.0021 


Experimental syst. 


±0.0067 


±0.0047 


±0.0056 


±0.0079 


Total hadronization 


±0.0035 


±0.0022 


±0.0023 


±0.0034 


a;^ = 0.5 


-0.0006 


+0.0001 


-0.0009 


<0.0001 


x^ = 2 


+0.0027 


+0.0016 


+0.0023 


+0.0016 


Total error 


+0.0089 
-0.0085 


+0.0066 
-0.0066 


+0.0070 
-0.0067 


+0.0089 
-0.0089 











7V<-' 


as {172 GeV) 


0.1046 


0.0990 


0.1075 


0.1028 


Fit range [logio(yc«t)] 


-2.75 0.75 


-2.5 0.75 


-2.75 1 


-2.75 1.25 


xVd-o.f. 


6.28/7 


2.53/6 


4.48/7 


0.61/6 


Statistical error 


±0.0046 


±0.0052 


±0.0037 


±0.0044 


Experimental syst. 


±0.0115 


±0.0107 


±0.0053 


±0.0062 


Total hadronization 


±0.0037 


±0.0016 


±0.0019 


±0.0026 


x^ = 0.5 


-0.0004 


+0.0002 


-0.0008 


-0.0001 


x^ = 2 


+0.0023 


+0.0015 


+0.0021 


+0.0015 


Total error 


+0.0131 
-0.0129 


+0.0120 
-0.0120 


+0.0070 
-0.0068 


+0.0082 
-0.0081 



Table 10: Results for as from fits of the In i^-matching predictions for the four observables 
, D^, and A^^ at i/s = 133, 161 and 172 GeV 
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a, (183 GcV) 


0.1076 


0.1019 


0.1100 


0.1048 


Fit raiig(- [loKio(/;,:a()] 


-2.75 -0.5 


-2.5 -U.3 


-2.75 -0.73 


-2.75 -0.75 


xVd-o.f. 


6.36/8 


9.8/7 


6.85/8 


5.07/8 


Statistical error 


±0.0030 


±0.0032 


±0.0021 


±0.0022 


Experimental syst. 


±0.0020 


±0.0052 


±0.0029 


±0.0046 


Total hadronization 


±0.0025 


±0.0014 


±0.002 


±0.0026 




-0.()0()5 


+0.()()02 


-0.0008 


-0.0001 


x^ = 2 


±0.0025 


+0.0017 


+0.0023 


+0.0015 


Total error 


+0.0050 
-0.0044 


+0.0063 
-0.0063 


+0.0047 
-0.0042 


+0.0059 
-0.0057 














^^(189 GeV) 


0.1109 


0.1062 


0.1106 


0.1079 


Fit range [logm (l/c«f)] 


-2.75 0.5 


-2.75 0.5 


-2.75 0.75 


-2.75 0.75 


\-/(l.().f. 


11.87/8 


5.12/8 


5.;-!(i/8 


(i.87/8 


Statistical error 


±0.0016 


±0.0015 


±0.0013 


±0.0015 


Experimental syst. 


±0.0029 


±0.0021 


±0.0043 


±0.0038 


Total hadronization 


±0.0022 


±0.0024 


±0.0017 


±0.0019 


Xn = 0.5 


-0.0005 


-0.0001 


-0.0008 


-0.0001 


x^ = 2 


+0.0028 


+0.0022 


+0.0023 


+0.0017 


Total error 


+0.0049 
-0.0041 


+0.0041 
-0.0035 


+0.0053 
-0.0049 


+0.0049 
-0.0045 



Table 11: Results for ag from fits of the In i?-matcliing predictions for the four observables 
, D^, and A^^ at ^/i = 183 and 189 GeV 







{y2y 


as{Mzo) 
xVd.o.f. 


0.1217 

2.77/7 


0.1162 

1.66/7 


Stat, error 


±0.0009 


±0.0006 


T()( al oxp. 


±0.0025 


±0.0021 


Total hadronization 


±0.0024 


±0.0037 


x^ = 0.5 


-0.0042 


-0.0031 


— 2 


+0.0056 


+0.0043 


Total error 


+0.0066 
-0.0054 


+0.0060 
-0.0052 



Table 12: Results for ^^(M^o) from fits of the 0{al) predictions for the energy evolution of 
(y2)^ and (ys)"' 



References 
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i?3'(2/™t = 0.08) 


R^{ycut = 0-01) 


R^ivcui = 0.01) 


X /d.o.j. 


0.1243 

3.Dd/7 


0.1127 

C 1 o /it 

5.18/7 


0.1137 

31.16/7 


fitted Xfjt 
stat. error 
xVd.o.f. 


U. 1 uou 

0.063 
±0.025 
2.98/6 


n 1 nQ7 

u. 1 uy i 

0.310 
±0.148 
5.03/6 


0.464 
±0.172 
30.70/6 


Stat. error 


±0.0003 


±0.0002 


±0.0003 


Total exp. 


±0.0024 


±0.0011 


±0.0014 


Total hadronization 


±0.0032 


±0.0082 


±0.0072 


x^ - 0.5 


-0.0061 


-0.0024 


-0.0014 


Xfj^ — 2 


+0.0074 


+0.0037 


+0.0030 


Total error 


+0.0085 
-0.0073 


+0.0091 
-0.0086 


+0.0079 
-0.0075 



Table 13: Results for ^^(M^o) from fits of the 0{a.l) predictions for the energy evolution of 
, and 
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